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Abstract 

A collective model is proposed to describe the chiral rotation and vibration and applied to a 
system with one /in/2 proton particle and one /in/2 neutron hole coupled to a triaxial rigid rotor. 
The collective Hamiltonian is constructed from the potential energy and mass parameter obtained in 
the tilted axis cranking approach. By diagonalizing the collective Hamiltonian with a box boundary 
condition, it is found that for the chiral rotation, the partner states become more degenerate with 
the increase of the cranking frequency, and for the chiral vibrations, their important roles for the 
collective excitation are revealed at the beginning of the chiral rotation region. 
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I. INTRODUCTION 



Since the chirality in nuclear physics was originally suggested by Frauendorf and Meng in 
1997 jl|, it has been one of the hot topics in nuclear physics. The corresponding experimental 
signals, the chiral doublet bands, were first observed in N = 75 isotones in 2001 and so 
far have been extensively reported in more than 30 nuclei. 

For a rotating nucleus with triaxiality, the three angular momenta, including collective 
rotor angular momentum along the intermediate axis, together with angular momenta of 
the valence particles (holes) along the nuclear short (long) axis, can construct the chiral 

n n 

geometry [1|, |2j. It has been demonstrated that the chiral doublet bands generally start 
from a left-right chiral vibration mode and evolve into a static chiral mode with the increase 
of spin 0,0]. 

Theoretically, chiral doublet bands were first predicted by the tilted axis cranking (TAC) 
approach and particle rotor model (PRM) [l]. Subsequently, numerous efforts have been 
devoted to the development of both the TAC methods 5j-[7| and PRM models % H, |9J. Both 
PRM and TAC have their own advantages and disadvantages. The PRM is a quantal model 
consisting of the collective rotation and the intrinsic single-particle motions, and the total 
angular momentum is a good quantum number. It describes a system in the laboratory 
reference frame where the spontaneously broken chiral symmetry in the intrinsic reference 
frame is restored. In the PRM, the energy splitting and quantum tunneling between the 
doublet bands can be obtained directly. However, rigid rotor with quadrupole deforma- 
tion parameters (3 and 7 has to be assumed at the very beginning in the practical PRM 
calculation. 

TAC approach is based on mean-field theory, and permits the calculation for the orien- 
tation of the density distribution relative to the angular momentum vector. However the 
TAC cannot describe the energy splitting and the quantum tunneling between the chiral 
doublet bands til. Up to now, the TAC method based on the Woods-Saxon or Nilsson po- 

n fiQ 

tential [5j, as well as the microscopic self-consistent Skyrme Hartree-Fock model p, LZj have 
been devoted to study the chirality. 

To describe the energy splitting between the chiral doublet bands, one has to go beyond 
the mean-field approximation. So far, this has been done in the framework of tilted axis 
cranking plus the random phase approach (RPA) model 0, [k]]. However, this model is 
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restricted in the description of the chiral vibration since the quantum tunneling in the chiral 
rotation is beyond the realm of RPA. Therefore, it is imperative to search a unified method 
for studying both chiral rotation and vibration in the framework of TAG. 

In analogy to Bohr Hamiltonian, constructing a collective Hamiltonian on the TAC so- 
lutions provides one of the ways out. Instead of the (3 and 7 degrees of freedom in Bohr 
Hamiltonian, a chiral degree of freedom should be introduced. Subsequently, the quantal 
tunneling in the region of chiral rotation can be described by considering the chiral fluctu- 
ations around mean-field minima besides the region of chiral vibration naturally described 
by the collective Hamiltonian. 

In this work the collective Hamiltonian for a system of one fan/2 proton particle and one 
fan/2 neutron hole coupled to a triaxial rigid rotor is constructed. The potential energy and 
mass parameters involved in the collective Hamiltonian are extracted from the TAC calcula- 
tions. By diagonalizing the collective Hamiltonian, the energy levels and wave functions are 
calculated and discussed in detail. The paper is organized as follows. In Sec. HIl based on 
the tilted axis cranking approach, the procedures for constructing and solving the collective 
Hamiltonian are introduced. The numerical details are presented in Sec. IIHI In Sec. IIVI 
the potential energy and mass parameters obtained from TAC approach are shown and the 
corresponding energy levels and wave functions obtained from the collective Hamiltonian 
are discussed in details. Finally, the summary is given in Sec. [V] 



II. THEORETICAL FRAMEWORK 



A. Tilted axis cranking model 



The detailed formalism for TAC can be found in Ref. 111. For a schematic discussion, 

n 

similar as Ref. pj, we consider a system of the fan/2 proton particle and the fan/2 neutron 
hole coupled to a triaxial rigid rotor. The cranking Hamiltonian reads 

fa' = fadcf - 3 ■ j, (1) 

3 = Jtt+Jv, (2) 

uj = (cj sin 9 cos <p, cj sin 9 sin ip, cu cos 0) , (3) 
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where the Hamiltonian of the deformed field is /idef = ^def + ^def with the single- j shell 
Hamiltonian 

ht? = ±\c{(jl - ^±12) cos 7 + ^=Ql sm 7 }. (4) 

The TAC solutions are obtained self-consistently by minimizing the total Routhian surface 

E\9 lV ) = {h')-\Y j J k ui (5) 

k=l 

with respect to the angles 9 and tp, where the moments of inertia for irrotational flow are 
adopted, i.e., 

J k = J oS in 2 ( 7 -y/c). (6) 

As discussed in Ref. there are three kinds of different solutions which can be dis- 
tinguished by the different values of 9 and tp: a) Principal axis cranking (PAC) solution: 
9 = 0,tt/2, tp = 0,±7r/2. b) Planar TAC solution: 9 ^ 0,±tt/2, tp = 0,±tt/2 or 9 = tt/2, 
tp ^ 0, ±7r/2. c) Aplanar TAC solution: 9 ^ 0, tt/2, tp ^ 0, ±7r/2, i.e., the chiral solution. 

In Ref. 12| , an orientation operator a was introduced to characterize the chiral degree of 
freedom. In the TAC approach for chiral solutions, the orientation operator is written as 

o" = (Jtt x j v ) ■R=\j ir \- \j v \ -\R\- sin 6> PN • sin 9 sin tp, (7) 

in which 6*pn denotes the angle between of j n and j u . This operator has opposite signs 
for the left-handed and the right-handed systems. Follow Ref. l|, we restrict the angle 9 
to < 9 < n/2, and by varying tp from —tt/2 to tt/2, we get tp as a dynamical variable 
describing a transition from the left-handed to the right-handed system. Thus for a given 
9, the angle tp can be used to characterize the chiral degree of freedom. We note that the 
calculations performed in Refs. and [h| for the potential energy as a function of 9 and 
tp for 134 Pr have shown a softer nature of the potential in the tp direction compare to the 9 
direction. 



B. Collective Hamiltonian 

Based on the TAC model, a collective Hamiltonian including the chiral degree of freedom 
could be constructed. The classical form of a collective Hamiltonian can be written in terms 
of tp as 

H con = T vih (p) + V(p), (8) 



where (f is the variable which characterizes the chirality, and the potential energy V(ip) 
could be extracted by minimizing the total Routhian Eq. fl5]) with respect to 9 for given tp. 
The vibrational kinetic energy of the chiral degree of freedom can be written as 

T vib = l -B^p\ (9) 

where B is the ma ss p arameter. The Hamiltonian flH]) is quantized according to the general 
Pauli prescription 13j, i.e., for a classical kinetic energy, 

T = \H B '^)m j , (io) 

the corresponding quantized form reads 

For the vibrational kinetic energy in Eq. ((9]), the mass parameter takes a 1 x 1 matrix form. 
Thus, the corresponding quantized form reads 

ft 2 d 1 d 

Hkin = ~2^W)^7W)^p' (12) 

Therefore, the quantized form of the collective Hamiltonian in Eq. dSJ) turns out to be 

h 2 d 1 d 
H coll = -2==^- — t==^- + V((p). (13) 

Note that the volume element in the present collective space is 



dr coU = / dip^B((p), (14) 



and the quantized Hamiltonian of Eq. f|T3|) is Hermitian with respect to the collective measure 
in Eq. (}TIj). 



C. Mass parameter 

The details for calculating the mass parameter can be found in fl^ . In this subsection, 
for completeness, a brief procedure is presented. Considering a variable ip as a function of 
time, the time dependent Schrodinger equation is written as 

h'(t)\m)=i^\m)- as) 



The general solution of this equation is 

m)) = Y. a ^ k(t) \ k )' 



(16) 



— / Ek(t')dt' and Ek(t) is the eigenvalue of the Hamiltonian h'(t). By 
h Jo 



where 4>k(t) = 

substituting Eq. ffl6|) into Eq. ffl5|) . we obtain 







oJ2a k (t)e^^(l\-^\k). 



(17) 



We assume that the wave function is the lowest eigenstate of the cranking Hamiltonian 
h' in Eq. dTJ at i = 0, i.e., a k = Sko at t — 0. Here, |0) means the lowest state of h! . If we 
further assume that a^o are small and the variable ip follows the relation of Cp = — Q 2 (p with 
Q being the vibrational frequency, the solution of the Eq. (fTTj) can be obtained as 

ih(Ei- E )(p + h 2 tt 2 (p d 



ai (E t - E ) 2 - h 2 tt 2 W dtp l ) 

where E\ and Eq are the eigen energies of the cranking Hamiltonian h! in Eq. (JXJ 
The energy of the system can be written as 



(18) 



E(t) = Eo + J2( E i~ E o)\ a i\ 



E + J2(Ei - E ) 
E + \B{^ 2 . 



[(E t - e ) 2 - h 2 n 



212 



(19) 



Neglecting the high-order terms (oc if 2 ), the mass parameter can be obtained as 



B(<p) = 2ft 2 ^ 



(E - E o r\(l\JL\0}\ 



212 



[(E t - e ) 2 - h 2 n j 

(^-^o)|(/|K^]|0)| 2 

[(Eh - e ) 2 - h 2 n 2 ] 2 ■ 



(20) 



By substituting Eq. ([T]) into Eq. (120]) and using the relation 

(/|[^dcf-^-/,^]|0) = ^(/|J|0), 



one obtains, 



B( V ) = 2h 2 J2 



(E t - E 



z^o 



[(E t - e ) 2 - h 2 n 2 ] 2 



(21) 



(22) 



Once the vibrational frequency Q is known, the mass parameter B(ip) can be then de- 
termined by Eq. (1221) . Normally, Q is determined as follows: a) for the chiral rotations, 
it is approximately taken as zero because the barrier penetration between the left-handed 
and right-handed states is low; b) for the chiral vibrations, the potential energy V(tp) can 
be approximated by the harmonic oscillator potential \K(p 2 with the corresponding spring 
coefficient K. The mass parameter is 

B=£. (23) 
Combining Eqs. (122]) and ( 123]) . the value of Q can be obtained for the chiral vibrations. 



D. Basis space 

After the mass parameter B(<p) is obtained from the TAC model, the collective Hamil- 
tonian (|T3|) is constructed. It is easy to find that the collective Hamiltonian Eq. (fT3j) keeps 
the parity conservation with respect to tp—t—ip. Therefore, the eigenstates of the collective 
Hamiltonian can be divided into two separate subspaces, i.e., the positive parity states and 
the negative parity states. 

For the positive parity, the basis states can be taken as 

, . . 12 cos(2n — l)tp „ 

and for the negative parity they are 

, . . /2sin2nw _ 

These basis states fulfill the box boundary condition, 

Vn(vr/2) = ^(-tt/2) = 0. (26) 

For the chosen basis, the normalization conditions with respect to the collective measure in 
Eq. (JUj) are satisfied. 

The wave function can be then expanded by the basis states as 

,. . *^ \2 cos(2n — l)ip /2sin2no9 . _ 

n=l v ' n=l 

where the expansion coefficients a n and b n {n > 1) are obtained by diagonalizing the collec- 
tive Hamiltonian Eq. ( 1T3|) . 
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III. NUMERICAL DETAILS 



In the following calculations, the symmetric particle-hole configuration 7r/i n / 2 <8> ^11/2 
is considered and 7 deformation is assumed as 7 = —30°. The coupling parameters C in 
the single-^ shell Hamiltonian are taken as C n = 0.25 MeV for the proton particle and 
C v = —0.25 MeV for the neutron hole, respectively. The moment of inertia is chosen as 
Jo = 40 h 2 /MeV. These parameters are the same as in Ref. [lj]. 



IV. RESULTS AND DISCUSSION 



A. Collective potential 

The potential energy surface in the rotating frame, i.e., the total Routhian E'(9,if) as 
a function of 9 and if, is shown in Fig. [T] at the frequencies tvjj = 0.1, 0.2, 0.3, 0.4 MeV. 
The present results are consistent with those in Ref. lj where the total Routhians surface 
calculations at the frequencies hu = 0.10, 0.30, 0.50 MeV have been presented. 

One can see that all the potential energy surfaces are symmetrical with the if = 0° 
line. This means that the two chiral configurations with ±\if\ for a given tilted angle 9 are 
identical. With the increasing frequency, the minima in the potential energy surfaces change 
from if = 0° to if ^ 0°. As discussed in Ref. [lj], this implies the rotating mode changes 
from planar to aplanar rotation. 

By minimizing the total Routhian E'(9, if) with 9 for given if, the potential energy V(if) 
in the collective Hamiltonian ffTB"]) is obtained and shown in Fig. [2j Again, the potential 
energy is symmetrical about if = 0° in correspondence with the results shown in Fig. [IJ For 



the frequency hu < 0.15 MeV, the potential V(if) is a harmonic osci 
only one minimum at if = 0°. This corresponds to the planar rotation 



lator type which has 
lj . For the frequency 



fkd > 0.20 MeV, the potential V(if) has two symmetrical minima. This corresponds to the 
aplanar rotation lj . Due to the appearance of the potential barrier which breaks the chiral 
symmetry, the stable chiral solutions are found in the body-fixed frame. The heights of 
barrier defined as AV = V(0) — Knin (in MeV) with Kn iT1 being the value of the potential 
at the minimum presented also in the figure. It is found that the potential barrier increases 
with the cranking frequency, e.g., from 1 keV at fko — 0.20 MeV to about 2 MeV at 
hcu = 0.50 MeV. 
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FIG. 1: (Color online) Total Routhian surface calculations for the h\\/2 proton particle and the 
h\\/2 neutron hole coupled to a triaxial rotor with 7 = —30° at the frequencies hw = 0.1, 0.2, 0.3, 
0.4 MeV. All energies are normalized with respect to the absolute minimum (star). The step is the 
energy difference between the contour lines. 



B. Mass parameter 

For chiral rotation with hw > 0.20 MeV, the chiral vibration frequency Q in Eq. ( 122]) is 
taken as Q = 0, which results in the mass parameter 

In Fig. [31 the mass parameter B(<p) as a function of tp for the chiral rotation cases is 
shown. It is seen that B(ip) is symmetric with respect to tp = 0° and increases dramatically 
when tp approaches to ±90°. In the interior part, B(tp) is increased remarkably with \tp\ for 
hu > 0.35 MeV, while its dependence on tp is weak for hw = 0.25 and 0.30 MeV. 
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FIG. 2: (Color online) The potential energy V(<p) as a function of tp extracted from the total 
Routhian surface calculations. The arrow labels the position of the potential minimum F m i n . The 
potential barriers defined as AV = V(0) — V m i n (in MeV) are also presented. 
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FIG. 3: (Color online) The mass parameter B{ip) as a function of tp for the chiral rotation cases 
obtained based on TAC. 



C. Energy levels 

The collective Hamiltonian can be constructed by Eq. (I13p with the potential energy V(<p) 
and mass parameter B(ip) obtained. The diagonalization of the collective Hamiltonian yields 
the energy levels and wave functions for each cranking frequency. As B(<p) is quite large at 
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ip ~ ±90°, a wall approximation is adopted. This corresponds to the choice of basis in Eq. 

In Fig. |U the six lowest energy levels, labeled as 1-6, obtained from the collective Hamil- 
tonian are shown together with the potential energy V((p). It is seen that with the increasing 
frequency, the three pairs of energy levels, i.e., levels 1 and 2, levels 3 and 4 as well as levels 
5 and 6, become close to each other. For example, the energy difference between levels 1 
and 2 decreases from 1.19 MeV at fiu = 0.25 MeV to 0.01 MeV at hu = 0.50 MeV. Since 
the potential barrier becomes higher and wider with the increase of cranking frequency, the 
tunneling penetration probability is more and more suppressed. Therefore, the two levels 
tend to be more degenerate. 
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FIG. 4: (Color online) The six lowest energy levels, labeled as 1-6, obtained from the collective 
Hamiltonian. The potential energy V(ip) is included as well. 



This can be more clearly seen in Fig. |5] where the energy difference AE between the 
lowest two levels of the collective Hamiltonian is shown as a function of the potential barrier 
height AV. It shows that with the increase of the cranking frequency, the potential barrier 
grows from 1 keV to about 2 MeV, while the energy differences AE drops correspondingly 
from 1.19 MeV to 0.01 MeV. 
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FIG. 5: (Color online) The energy difference AE between the lowest two levels 1 and 2 of the 
collective Hamiltonian as a function of the potential barrier height AV . 
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FIG. 6: (Color online) Wave functions ip(^p) for the lowest two levels 1 and 2 obtained from collective 
Hamiltonian. 



D. Wave function 

The wave functions ip(<p) for the lowest two levels are shown in Fig. [6] and the corre- 
sponding probability distributions determined as |i/7(</?)| 2 are shown in Fig. [TJ It is found 
that the wave functions are symmetric for level 1 and antisymmetric for level 2 with respect 
to ip — > —ip transformation. Thus the chiral symmetry broken in the aplanar TAC solutions 
is restored. It is also shown that the state with symmetric wave function is favored in energy 
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FIG. 7: (Color online) Probability distributions for the lowest two levels 1 and 2 calculated by 



compared with the corresponding state with antisymmetric wave function. 

For the frequency hw = 0.25 and 0.30 MeV, the wave function of the level 1 remains 
almost constant except at the border ip = ±90° which indicates that the fluctuation of the 
total angular momentum is large. The wave function of the level 2 shows peaks close to 
(f = ±90°, which suggests that the chiral vibration plays an important role for the collective 
excitation at the beginning of the chiral rotation region. When the cranking frequency 
increases, the wave function of level 1 tends to show similar pattern. This implies the 
appearance of the nearly identical physical properties of these two levels and a good static 
chirality appears. 

It is well known that the TAC results are given in the body-fixed frame and the two chiral 
partner solutions are identical. It cannot describe tunneling between the left- and right- 
oriented solutions. The chiral symmetry is thus broken and only one rotational band can 
be obtained. However, in the present formulation of the collective Hamiltonian, the results 
are given in the laboratory reference frame, which could describe the quantum tunneling 
between the two chiral solutions. Moreover, the chiral symmetry is restored in the present 
framework by considering the contributions from all the <p directions in the potential energy 
surface. 
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E. Comparison with exact solutions 
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FIG. 8: (Color online) The energy spectra of the doublet bands obtained from the collective 
Hamiltonian in comparison with the exact solutions by the PRM. Inset: The spin I(uS) obtained 
from TAC in comparison with PRM. A similar I (to) plot has been given in Ref. [lj]. 



The exact solutions for the system discussed here can be obtained by the particle rotor 
model. In Ref. jl], the relation between the angular momentum / and the rotational fre- 
quency hjj as well as the intra band transition probabilities obtained from the TAC has been 
compared with the PRM results. It has been shown that the TAC solutions could reproduce 
the results of the yrast levels of the PRM calculation jl]]. Here, in the inset of Fig. 8, the 
spin obtained from TAC in comparison with PRM is shown once more, where the good 
agreement can be clearly seen. 

In Fig. |Hl the energy spectra of the doublet bands obtained from the collective Hamilto- 
nian based on TAC are compared with the PRM results. One can see that apart from the 
agreement of collective Hamiltonian and PRM results for the yrast band, the partner band 
of PRM can also be reasonably reproduced by the collective Hamiltonian. 

For the collective Hamiltonian results, the energy differences between the doublet bands 
become smaller with the increase of the cranking frequency. For the PRM results, however, 
the doublet bands become closer up to hu ~ 0.35 MeV and the energy differences between 
the doublet bands continue to increase for the higher cranking frequency. As demonstrated 
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14, 



both in PRM |4j, |l5j, [16| and TAC+RPA [10[ investigations, the doublet bands will attain a 



second chiral vibration character, which is not taken into account by the present collective 
Hamiltonian investigation. 

The present comparison for the energy spectra is made with respect to rotational fre- 
quency hut, rather than the observable of spin I. In the TAC approach, the angular momen- 
tum is not a good quantum number, and the diagonalization of the collective Hamiltonian 
is carried out for a given rotational frequency. Therefore, it is appropriate to compare the 
energy spectra with respect to the rotational frequency with the results of PRM. It can be 
seen from the inset of Fig. |H]that, for a certain spin, the rotational frequencies of the doublet 
bands 1 and 2 calculated by PRM may be quite different (for example, the differences are 
respectively 0.06, 0.0, 0.02, and 0.04 MeV at I = 12,16,20,24 h). Here the fluctuation of 
the potential energy with 9 is neglected and only the if is treated as a dynamical variable 
describing the transition from the left-handed to the right-handed system. The success of 
the collective Hamiltonian here guarantees its application for realistic TAC calculations. 



V. SUMMARY 



In summary, a collective model which is able to describe the chiral rotation and vibration 
is proposed and applied to a system of one hn/ 2 proton particle and one /in/2 neutron hole 
coupled to a triaxial rigid rotor. In this framework, it goes beyond the mean-field approxi- 
mation, includes the quantum fluctuation in the chiral degree of freedom, and restores the 
chiral symmetry. Based on the tilted axis cranking approach, both the potential energy and 
mass parameter as functions of ip are obtained and included in the collective Hamiltonian. 
Diagonalizing the collective Hamiltonian with a box boundary condition, the energies and 
the wave functions of the collective states corresponding to the motion along the chiral de- 
gree of freedom are obtained. It is found that for chiral rotation, the partner states become 
more degenerate with the increase of the cranking frequency, and for the chiral vibrations, 
their important roles for the collective excitation are revealed at the beginning of the chiral 
rotation region. 
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